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TOWARD  IN-SITU  BUILDING 
R-VALUE  MEASUREMENT 


Stephen  N.  Flanders  and  Stephen  J.  Marshall 


INTRODUCTION 

Measurement  of  building  envelope  thermal  per¬ 
formance  reveals  the  potential  for  investment  in 
improved  insulation,  in  both  existing  and  new  con¬ 
struction.  We  are  developing  a  technique  for  such 
measurement.  It  combines  localized  thermal  meas¬ 
urements,  using  heat  flow  and  temperature  sensors, 
with  infrared  (IR)  thermal  mapping  to  assess  the 
R-values  of  large  areas.  In  this  paper  we  summarize 
our  efforts  to  date,  which  encompass  measurements 
of  wood  frame  and  masonry  buildings.  The  technique 
that  we  and  other  investigators  have  employed  is  rea¬ 
sonable,  but  requires  further  exploration  so  that  all 
factors  affecting  accuracy  are  better  controlled. 

In  brief,  an  IR  scanner  first  surveys  the  area  to 
determine  the  apparent  warmest  and  coolest  loca¬ 
tions  on  the  building  surface.  Contact  thermal  sen¬ 
sors  at  these  locations  then  provide  heat  flow  and 
temperature  data  that  determine  the  areas  of  highest 
and  lowest  thermal  resistance  (R-value).  A  thermo¬ 
graphic  map  that  depicts  all  the  gradations  of  tem¬ 
perature  between  the  instrumented  locations  is  pre¬ 
pared,  forming  the  basis  for  interpolating  R-values 
for  each  location.  We  outline  this  technique  in 
Flanders  and  Marshall  (1981). 

The  investigation  reported  here  encompassed 
field  measurements  of  masonry  and  frame  construc¬ 
tions.  The  masonry  construction  was  a  cavity  wail 
of  brick  and  concrete  blocks.  The  wood  frame  con¬ 
struction  incorporated  full-thickness  insulation  be¬ 
tween  89-mm-deep  2x4  studs  and  raftets.  Computer 
simulations  of  overall  wall  behavior  and  of  heat  flow 


around  the  sensors  provided  an  independent  check 
on  the  field  results.  Laboratory  investigations  dem¬ 
onstrated  the  different  effects  of  sensor  configuration 
and  calibration  methods. 

Other  investigators  (for  example,  Treado  1980) 
have  achieved  good  agreement  between  measured 
and  calculated  values  for  thermal  performance  of 
building  envelopes.  On  the  other  hand,  our  work 
with  masonry  construction  (Flanders  and  Marshall 
1982)  resulted  in  measured  R-values  that  were  57% 
of  those  expected  for  locations  containing  no  insula¬ 
tion.  The  large  discrepancies  cast  doubt  on  how  well 
the  results  reflected  the  actual  performance  of  masonry 
construction.  Our  measurement  of  four  walls  and 
two  roof  surfaces  of  a  frame  building,  reported  here 
for  the  first  time,  gave  us  measured  R-values  only  13% 
less  than  the  theoretical  value.  These  results  suggest 
difficulties  with  measuring  masonry  construction  and 
that  the  technique  may  be  reasonably  accurate  for 
frame  construction. 

R-VALUE  MEASUREMENTS  ON  A 
FRAME  BUILDING 

A  4.2-  by  7  J-m  frame  building  with  a  1 2.5°  pitched 
roof  was  instrumented  with  heat  flow  sensors  (HFS’s) 
and  thermocouples  on  each  of  four  walls  and  the  two 
halves  of  the  low-pitch  gable  roof  (see  Fig.  I ).  Simul¬ 
taneous  recording  of  data  from  these  locations  dem¬ 
onstrated  the  effect  of  orientation  of  construction 
on  the  measurement  (the  building  axes  are  rotated 
18s  clockwise  of  the  true  compass  axes).  Figure  2 
shows  the  typical  roof  and  wall  construction,  which 


Figure  1.  The  frame  building  tested  (top).  The  site  is  sheltered  between  two 
ridges  against  prevailing  winds,  which  are  light.  Sensor  locations  on  the  building 
(bottom). 


incorporated  89  mm  of  insulation  throughout.  Each 
of  the  six  sensor  locations  was  equidistant  between 
framing  members. 

Each  sensor  location  employed  thermocouples  on 
indoor  and  outdoor  surfaces  and  an  HFS  on  the  in¬ 
door  surface.  Each  outdoor  wall  thermocouple  had 
a  1 2  by  1 2  mm  gable  roof  projecting  over  it  for  shade 
against  direct  sunlight.  A  strip  of  vinyl  covered  the 
thermocouple  on  the  roof. 

We  used  two  types  of  HFS’s,  one  with  and  one 


without  a  built-in  guard  surrounding  the  sensing 
area.  The  guarded  sensor  incorporated  a  thermopile 
measurement  area  of  25  by  13  mm  within  an  over¬ 
all  area  of  73  by  1 1 1  mm  and  it  was  3.2  mm  thick. 
The  material  was  an  unidentified  laminated  plastic 
resin.  The  unguarded  sensor  had  a  shiny  black-coated 
layer  of  aluminum  bonded  to  a  dull  black  plastic 
sheet  containing  the  thermopile.  Its  overall  dimen¬ 
sions  were  13  by  51  by  2.5  mm. 

The  HFS’s  were  calibrated  in  a  Rapid-k  heat  flow 
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Figure  2.  Typical  construction  of  the  frame  building,  incorporating  2  *4s  spaced  40.6  cm 
on  center  in  both  the  roof  and  wall. 


meter  thermal  conductivity  testing  machine,  which 
maintains  set  temperatures  on  the  upper  and  lower 
faces  of  a  0.3-m  specimen  and  reads  heat  flow  from 
a  100-mm  square  region  in  the  center.  The  test 
specimen  had  three  0  J-m  square  layers  of  material 
with  a  total  resistance  of  0.8S  K-m2/W;  the  sensor 
was  embedded  in  the  middle  layer.  The  top  layer 
was  6  mm  of  neoprene  foam  and  the  bottom  layer 
25  mm  of  rigid  fiberglass  insulation.  Between  these 
was  a  3 -mm  sheet  of  Mi  carta  (a  paper-reinforced 
phenolic  resin)  with  a  cutout  in  the  center  to  accom¬ 
modate  the  HFS.  The  calibration  temperatures  were 
chosen  to  approximate  the  temperature  that  the  HFS 
would  function  at  when  in  use. 

The  HFS’s  were  taped  to  the  interior  surface  with 
a  layer  of  petroleum  jelly  between  the  surface  and 
sensor  to  provide  thermal  contact.  The  tape  surrounded 
the  edge  of  die  sensor,  overlapping  it  by  about  5  mm. 
The  indoor  thermocouple  wu  taped  adjacent  to  die 
HFS.  A  sheet  of  white  paper  covered  the  indoor  sen¬ 
sors  at  each  location  to  damp  fluctuations  in  heat 
flow  readings  due  to  convection  and  to  mask  the 
sensor's  absorptivity  to  Mend  with  its  surroundings. 

A  thermographic  survey  at  a  Af  (difference  be¬ 
tween  indoor  and  outdoor  surface  t  tiperature)  of 
about  20*C  revealed  terms!  mattes  in  the 


Stud  bays,  despite  the  presence  of  surface  water  stains. 

In  particular,  there  was  no  evidence  of  air  leakage  or 
convection  in  the  plane  of  the  construction. 

We  began  recording  data  on  19  March  1981  at 
30-minute  intervals  and  stopped  after  750  hours. 

During  this  time  the  outdoor  dry  bulb  temperatures 
averaged  maximums  of  13°C  each  day  and  minimums 
of  -2.4°C  with  a  3°C  warming  trend.  The  wind  aver¬ 
aged  over  each  day  1 .5  m/s  from  die  west  with  a  maxi¬ 
mum  average  speed  for  one  day  of  3.5  m/s.  Twelve  days 
with  precipitation  were  2  ft  days  apart  on  the  average. 

Hie  R-value  was  calculated  by  summing  all  AT 
data  and  dividing  by  the  sum  of  all  heat  flux  data 


R  -  ZATfLQ 
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where  Q  is  the  heat  flux.  Equation  1  is  die  ratio  of 
average  values  for  AT  and  Q  which,  as  the  constant 
terms  in  a  Fourier  series  of  time-variant  data,  approx 
imste  the  definition  of  It  under  steady -state  condi¬ 
tions  whan  the  measurement  lasts  long  enough  for 
the  equation  to  converge  on  a  steady  value.  Both 
Tteado  (1980)  and  lohaonesson  (1979)  advocate 
the  use  of  eq  1  or  its  equivalent  in  integrals  derived 
ia  Appendix  A.  Figure  3  shows  a  typical  data  plot 
with  the  resulting  plot  of  It,  according  to  eq  1 . 


R  Value 
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0 
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Table  1 .  R -values  (K-m2/W)  calculated  according 
to  eq  1  for  7S0  hr  of  data  obtained  from  frame 
building.  Expected  R  value  *  2.29  K-m2/W  for 
all  measurement  sites. 


Wall  Roof 

Orientation  Orientation 

Measure _ N  E  S  W _ JV  t 


Rvalues  2.1  1.8  2.0  2.0  2.2  1.9 

Percent  difference  -9  -23  -13-13  -9  -19 

from  expected 


Table  2.  Consistency  of  R  value  readings  taken 
in  150-hr  blocks  expressed  as  a  percentage  of  the 
cumulative  750-hr  value.  Values  greater  than 
100%  were  more  than  the  750-hr  value. 

Watt  Roof 

Period  Orientation  Orientation 


(hr) 

N 

£ 

S 

W 

N 

S 

0-1  SO 

99 

S3 

99 

99 

89 

84 

150-300 

103 

124 

103 

106 

103 

104 

300-4SO 

97 

99 

96 

101 

108 

98 

450-000 

101 

102 

100 

102 

109 

106 

600-750 

99 

108 

99 

97 

97 

106 

Avnao«% 

100 

103 

99 

101 

101 

100 

Std.  dev. % 

2 

14 

2 

3 

8 

9 

Table  I  summarizes  the  results  for  walls  and  roofs 
with  the  same  expected  R-vaiue  of  2.29  K  m2/W. 

Ail  R-values  calculated  from  the  data  averaged  1 3% 
less  than  the  expected  value  based  on  the  thermal 
resistances  of  the  constituent  materials  listed  in 
ASHRAE  (1977). 

Dividing  the  data  into  five  groups  of  readings  1 50 
hr  in  duration  caused  about  the  same  percentage 
difference  from  the  expected  value,  but  the  scatter 
of  calculated  R-values  had  a  standard  deviation  that 
ranged  from  ±  2  to  1 1  percentage  points  around  the 
percent  difference  shown  in  Table  1. 

The  data  for  each  150-hr  block  were,  on  the 
whole,  repeatable.  The  R-value  obtained  during  any 
one  period  for  most  locations  was  within  a  standard 
deviation  of  10%  of  the  long-term  750  hr  value,  as 
Table  2  shows. 

The  deviations  of  calculated  R-values  from  ex¬ 
pected  shown  in  Table  1  probably  reflect  the  accur¬ 
acy  of  the  measurement,  based  on  our  thermography, 
which  detected  no  air  movement  at  Af=  26.1  °C, 
and  based  on  an  independent  check  of  the  R-value 
of  the  frame  construction.  A  0  J-m  square  specimen 
of  fiberglass  wool  insulation  was  obtained  from  di¬ 
rectly  behind  the  south  wall  sensor  location  (Fig.  1). 
The  specimen  was  placed  in  a  Rapid-k  thermal  con¬ 
ductivity  testing  machine.  The  R-value  obtained 
from  this  steady-state  laboratory  test  was  only  4% 
less  than  the  value  listed  in  ASHRAE  (1977). 
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Since  the  plywood  skins  of  the  construction  con¬ 
stitute  such  a  small  portion  of  the  overall  thermal 
resistance,  they  could  not  be  responsible  for  the  dif¬ 
ference  between  the  measured  and  expected  R-values 
in  the  Held  tests.  The  sensors  were  located  where 
heat  flow  should  be  perpendicular  to  the  plane  of 
the  wall.  Therefore,  we  must  ascribe  the  discrepancies 
to  the  measurement  process.  The  problem  is  to  iden¬ 
tify  which  factors  are  important  sources  of  error  and 
how  to  correct  them.  Although  Table  1  indicates 
that  the  measurement  technique  is  somewhat  biased, 
Table  2  demonstrates  that  the  measurements  are  re¬ 
peatable. 

FACTORS  AFFECTING  ACCURACY 
OF  MEASUREMENT 

There  are  several  possible  sources  of  inaccuracy  in 
the  measurement  and  calculation  of  R-values,  including: 

•  Insufficient  combination  of  AT  and  measure¬ 
ment  duration  for  a  converged  calculation, 
according  to  eq  1 . 


•  Multidimensional  heat  flow  around  the  HFS, 
and  lateral  heat  flow  within  die  construction 
due  to  conduction  or  air  movement . 

•  Inadequate  thermal  contact  of  the  HFS  or 
thermocouple  with  the  construction. 

•  Unrepresentative  exposure  of  sensors  to  sources 
of  infrared  radiation. 

The  effects  of  time  and  AT 

The  frame  building  was  measured  for  more  than 
750  hr.  With  an  average  AT  of  about  1 1®C,  this 
would  give  ample  time  for  the  proportion  of  unre¬ 
solved  transient  heat  flow  responses  ( Q )  to  varying 
AT  inputs  that  have  not  resulted  in  a  heat  flow  re¬ 
sponse  on  the  indoor  surface  to  be  very  small.  There¬ 
fore,  we  took  the  cumulative  R  values  after  750  hr 
to  be  the  reference  values  for  comparing  with  the 
R-values  obtained  during  each  150-hr  interval.  Not 
only  did  the  value  of  R  obtained  during  each  interval 
vary,  as  mentioned  above,  but  the  period  necessary 
for  it  to  stabilize  within  10%  of  the  750-hr  value 
varied,  from  less  than  10  up  to  76  hr  for  walls.  Roofs 
required  longer  times,  as  Figure  4  demonstrates. 


Figure  4.  Duration  of  mamrement  necessary  to  obtain  an  R -value  with¬ 
in  10%  of  the  value  obtained  after  750  hr  as  a  function  of  averaged  AT. 
Bach  set  of  data  for  an  average  AT  commences  at  the  beginning  of  the 
first  of  five  consecutive  150-hr  periods. 
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Time  to  stabilize  within  1 5%  of  the  750-hr  R-value 
was  not  a  strong  function  of  AT  within  the  range  of 
AT  values  we  observed. 

Our  work  with  masonry  walls  (Flanders  and  Mar¬ 
shall  1982)  encompassed  a  much  wider  range  of  AT 
values  and  demonstrated  that  AT  is  the  most  power¬ 
ful  variable  affecting  the  speed  of  convergence  at  an 
R-value  (calculated  with  eq  1).  In  that  paper  we  re¬ 
ported  that  extending  the  measurement  of  R  for  a 
masonry  wall  to  473  hr  changed  the  calculated  result 
by  only  8%  of  the  value  obtained  after  1 16  hr. 

Multidimensional  heat  flow 

Multidimensional  heat  flow  around  an  HFS  can 
result  both  from  the  nature  -/f  the  construction  and 
from  how  the  sensor  itself  affects  the  area  subject  to 
measurement.  The  masonry  wall  we  reported  on 
(Flanders  and  Marshall  1982)  had  a  variety  of  paths 
that  would  keep  heat  from  flowing  perpendicularly 
through  the  construction  (shown  in  Fig.  5).  Masonry 
ties  concentrated  the  heat  flow  between  the  brick 
and  concrete  block  sides  of  the  construction.  None¬ 
theless,  applying  the  HFS  at  different  places  along 
the  face  of  the  concrete  block  resulted  in  at  most  an 
8%  variation. 

The  effect  of  the  HFS  on  the  surface  being  measured 
must  be  defined.  Schwerdtfeger  (1970)  illustrates  the 
importance  of  the  comparative  thermal  conductivities 
of  an  HFS  and  the  homogeneous  medium  surrounding 
it  (he  was  primarily  interested  in  heat  flow  through 
soils).  Given  the  geometry  of  an  HFS,  he  demonstrated 


that  heat  flow  through  the  sensor  can  be  as  much  as 
1 .8  or  as  little  as  0.2  that  through  the  surrounding 
material.  This  varies  as  the  ratio  of  the  two  thermal 
conductivities,  *M0Ior/*ni.,eri*|.  ran«cs  from  01  to 
10. 

Heat  flow  through  a  budding  surface  represents  a 
different  situation  from  Schwerdtfeger’s.  The  HFS 
shape  and  conductivity,  the  conductivity  of  the 
building  surface,  and  the  effect  of  air  convecting  on 
the  back  side  of  the  sensor  become  factors  affecting 
whethet  heat  flow  is  perpendicular  through  the  sur¬ 
face.  To  assess  this  effect,  we  conducted  experi¬ 
ments  in  a  large-scale  calibrated  hot  box  and  simu¬ 
lated  heat  flow  with  a  finite-difference  computer 
program. 

The  calibrated  hot  box-type  thermal  testing  box 
consisted  of  two  constant  temperature  chambers, 
one  refrigerated  and  the  other  heated,  on  either 
side  of  a  3.7-  by  2.4-m  test  wall.  The  wall  was  made 
of  extruded  polystyrene  foam,  50.8  mm  thick.  Four 
sensors  (two  guarded  and  two  unguarded  HFS's  as 
described  earlier)  were  installed  midway  up  the  wall. 
A  0.3-m  square  by  3-mm  thick  Micarta  guard  had  a 
cutout  in  the  center  to  accommodate  the  two  sizes 
of  HFS.  The  thermal  resistance  of  the  Micarta  and 
the  two  types  of  HFS  were  unknown,  but  are  be¬ 
lieved  to  be  similar.  Table  3  demonstrates  the  ef¬ 
fect  of  the  guard. 

The  net  effect  of  having  the  guard  for  either  sen¬ 
sor  type  was  a  5%  reduction  in  heat  flow  through 
the  HFS.  Since  the  heat  flow  through  the  sensor 


Gap  Between 
Brick  and  CMU 


UF  Foam 
(Discontinuous) 

Figures.  Configuration  of  masonry  well studied.  Webs  in  concrete  Mocks  end  masonry 
lift  divert  heat  flow  from  a  path  perpendicular  to  the  construction. 


Table  3.  Ratio  of  heat  flow  measured  with  the 
HFS  and  that  expected  from  the  mater¬ 

ial  properties  and  dimensions  (Gexp);  italics  de¬ 
note  a  reading  with  a  guard. 


Small  HFS  Large  HFS 

Senior  no.  Senior  no. 


Measure 

i 

2 

3 

4 

^meas^exp 

1.10 

1.04 

0.94 

1.02 

1.02 

1.02 

0.99 

0.99 

1.10 

0.91 

0.99 

0.99 

1.16 

0.99 

0.97 

0.9S 

Average,  no  guard 

1.10 

0.98 

1.00 

1.02 

Guard/a  ve.  no  guard 

0.93 

0.97 

0.96 

0.90 

without  the  guard  was  very  close  to  that  expected, 
the  guard  constituted  no  benefit  for  that  combina¬ 
tion  of  HFS  and  building  material. 

Another  test  was  designed  to  test  the  extremes 
to  which  an  HFS  might  distort  the  results  when  the 
surface  material  was  much  more  highly  conductive 
than  the  HFS.  This  test  employed  a  0.3-m  square 
specimen  consisting  of  a  6.4-mm  thick  aluminum 
plate  and  a  25.4 -mm  thick  extruded  polystyrene 
board.  These  were  placed  in  the  Rapid-k  thermal 
testing  machine  without  closing  the  internal  platform 
that  normally  presses  the  specimen  against  the  upper 
and  lower  surfaces  of  the  test  chamber.  Instead,  an 
HFS  placed  on  top  of  the  specimen  and  covered  with 
paper  was  adjacent  to  an  air  space. 

With  the  polystyrene  as  the  top  material  of  the 
specimen,  the  HFS  affixed  to  it  indicated  9%  more 
heat  flow  than  expected  for  the  test  setup.  With 
the  aluminum  as  the  top  material,  two  runs  gave 
readings  that  were  69%  higher  than  expected  from 
the  HFS.  This  test  requires  further  refinement  be¬ 
cause  of  the  high  lateral  conduction  of  the  aluminum, 
but  it  demonstrates  the  order  of  magnitude  of  dis¬ 
crepancy  possible  when  the  HFS  monitors  a  highly 
conductive  material. 

To  test  the  extent  to  which  an  HFS  disturbs  heat 
flow  on  a  wooden  surface  backed  with  fiberglass,  a 
computer  simulation  employed  the  finite-difference 
approach  to  solving  the  steady-state  condition.  A 
grid,  consisting  of  1 .27-mm  square  nodes,  89  nodes 
across  and  18  nodes  thick,  simulated  half  of  the 
long  dimension  of  an  HFS  in  two  dimensions.  Semi¬ 
infinite  nodes  surrounded  two  sides  of  this  region. 
Along  the  line  of  symmetry  there  were  flux  nodes 
with  no  lateral  flux.  The  side  with  the  HFS  was 
represented  with  convective  nodes  adjacent  to  the 
exposed  surfaces  ofwaB  and  sensor.  Nodes  at  top 


and  bottom  of  the  grid  were  a  constant  temperature 
at  a  semi-infinite  distance  away. 

The  simulation  indicated  that  heat  flow  at  the 
center  of  the  sensor,  which  represented  the  larger 
of  the  two  sensors  we  employed,  might  be  93%  of 
the  flow  where  there  was  no  HFS.  A  similar  analysis 
with  the  HFS  on  a  concrete  block  resulted  in  a  read¬ 
ing  of  99%  of  the  undisturbed  value. 

We  concluded  that  the  guard  does  not  contribute 
a  significant  benefit  for  thermal  measurements. 

Thermal  contact  and  convection 

Johannesson  (1979)  reported  that  a  1-mm  or 
smaller  air  gap  between  the  sensor  and  the  building 
surface  has  a  minimal  effect  on  sensor  accuracy.  We 
tested  a  commercial  heat-conductive  compound  and 
a  petroleum  jelly  to  prevent  convection  behind  the 
HFS  and  to  improve  thermal  contact.  We  also  cov¬ 
ered  the  indoor  sensors  with  paper  to  minimize 
fluctuations  in  the  readings.  Both  tests  employed 
the  large-scale  calibrated  hot  box  described  earlier. 

Hie  heat-conductive  compounds  had  no  signifi¬ 
cant  effect  on  sensor  readings.  The  HFS  readings 
with  petroleum  jelly  underneath  the  sensors  were 
99%  of  those  without  the  jelly ,  with  a  standard  de¬ 
viation  of  2%.  The  commercial  compound  caused 
readings  102%  of  the  base  case  with  a  standard  de¬ 
viation  of  5%. 

The  combination  of  paper  and  petroleum  jelly 
did  improve  the  constancy  of  the  output  from  die 
HFS’s.  Whereas  the  base  case,  demonstrated  in 
Table  3,  resulted  in  a  standard  deviation  of  2  to  5% 
of  the  heat  flux  reading,  the  addition  of  a  paper 
cover  and  petroleum  jelly  between  sensor  and  wall 
diminished  the  variability  to  between  2  and  3.5%,  an 
average  reduction  of  the  basic  standard  deviation  by 
about  27%. 

Exposure  of  sensors  to  sources  of  R  radiation 

The  R-values  reported  in  Table  1  indicate  that 
orientation  of  the  surface  with  respect  to  the  sun 
affects  the  measurement  accuracy.  The  sensors  on 
the  north-facing  wall  and  roof  gave  readings  only  6 
and  4%  less  than  the  theoretical  expected  value. 
However,  with  a  12.5°  roof  pitch  and  a  52s  average 
angle  between  the  sun’s  ecliptic  and  die  horizon, 
the  northerly  roof  surface  received  more  solar  radia¬ 
tion  than  did  the  northerly  wad.  The  readings  from 
the  other  orientations  were  12  to  23%  less  than  the 
expected,  although  the  south-facing  surfaces  did  not 
necessarily  have  die  greatest  variation. 

If  die  outdoor  temperature  sensors  were  reading 
higher  than  the  temperature  of  the  surface  to  which 
they  ware  attached,  then  die  AT  between  indoor  and 


outdoor  thermocouples  would  be  less  than  the  AT  be¬ 
tween  indoor  and  outdoor  surfaces.  This  would  cause 
the  R-value  calculation  from  eq  1  to  be  less  than  the 
actual  value. 

The  indoor  sensors  are  similarly  vulnerable  to 
inaccuracy  from  absorbing  thermal  radiation  differ¬ 
ently  than  the  surface  they  are  supposed  to  monitor. 

A  nominal  solution  to  this  problem  is  to  select  a 
cover  that  has  the  same  absorptivity  as  the  interior 
finish  of  the  construction.  Since  the  absorptivity  of 
most  materials  at  near  room  temperature  radiation 
is  about  0.9,  a  sheet  of  paper  should  suffice. 

A  few  cautionary  notes,  however:  the  absorptivity 
of  some  pastel-colored  paints  may  be  as  low  as  0.3 
(Marshall  1981).  Likewise,  the  absorptivity  of  white 
paper  can  range  from  0.7  to  0.9.  Furthermore,  even 
if  the  cover  does  not  touch  the  sensors,  some  heat 
from  the  cover  may  be  transferred  to  the  sensors  via 
radiation,  thereby  diminishing  the  value  of  the  cover 
as  a  way  to  eliminate  variation  in  absorptance. 

Measurement  improvements  needed 

Ftobably  the  most  significant  source  of  error  in 
the  measurements  of  the  wood-frame  construction 
was  the  outdoor  sensors,  which  registered  higher  than 
representative  temperatures.  Further  work  will  de¬ 
termine  how  severe  this  effect  is  and  how  to  dimin¬ 
ish  it. 

The  exterior  sensors  also  probably  contributed  to 
the  discrepancies  reported  in  our  ASHRAE  paper 
(Flanders  and  Marshall  1982)  on  measuring  masonry 
wall  thermal  conductivity.  However,  lateral  heat 
flow  through  the  masonry  construction  may  be  a 
bigger  cause  for  error  than  we  surmised.  Furthermore, 
convection  within  the  masonry  construction  may  have 
worsened  its  thermal  performance  more  than  expected. 


CONCLUSION 

One  reason  for  addressing  the  accuracy  of  thermal 
measurements  is  to  gain  an  understanding  of  the 
variables  that  will  permit  the  development  of  re¬ 
commended  practices  for  the  use  of  heat  flow  sensors 
and  thermocouples  for  characterizing  the  thermal 
performance  of  building  envelopes.  Currently,  an 
investigator  cannot  rely  on  the  literature  for  guidance 
in  assessing  the  limitations  of  such  thermal  measure¬ 
ments  or  in  obtaining  a  desired  level  of  accuracy. 

Such  a  recommended  practice  will  address  the 
use  of  sensors,  sensor  placement,  acceptable  thermal 
conditions,  and  measurement  duration,  among  o flier 
considerations.  Some  of  the  most  important  recom¬ 
mendations  are  listed  below. 


Sensors 

•  Calibrate  heat  flow  sensors  under  conditions 
that  correspond  to  those  of  the  test;  e.g.  for  a 
frame  wall  test,  calibration  in  a  vertical  orien¬ 
tation  with  gypsum  board  on  one  side  and  air 
on  the  other. 

•  Mask  all  sensors  with  a  material  of  the  same 
absorptivity. 

Sensor  location 

o  Avoid  known  areas  of  lateral  heat  flow  from 
conduction  and  air  flow 
o  Employ  thermography,  where  possible,  to 
identify  thermal  anomalies. 

Thermal  conditions 

•  Identify  tradeoffs  between  AT  and  duration 
of  measurement. 

Duration  of  measurement 

•  Determine  when  sufficient  data  have  been  ob¬ 
tained. 

For  now,  it  is  the  accuracy,  not  the  repeatability, 
of  the  technique  that  needs  improvement. 
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APPENDIX  A.  DERIVATION  OF  EQUATION  I 

There  are  two  avenues  toward  understanding  the  justification  for  eq  i ,  one  intuitive  and  the 
other  mathematical.  First,  we  will  consider  the  intuitive  point  of  view. 

R-value  is  actually  defined  only  for  steady -state  conditions  for  temperature  and  heat  flux: 

R  =  A  TIQ  (A 


where  AT  is  the  difference  in  temperature  between  the  warm  and  cool  sides  of  the  test  specimen 
and  £is  the  heat  flow  measured  coming  through  the  specimen.  This  simply  states  that  every  unit 
of  temperature  drop  across  the  construction  causes  a  proportionate  measure  of  heat  flow  through 
it,  reflecting  a  constant  property  of  the  constituent  materials.  This  concept  remains  true  for  the 
sum  of  all  time-varying  AT’s  across  a  construction  when  compared  with  the  sum  of  all  time-varying 
Q'%  passing  through  it,  both  taken  over  a  long  period.  The  shorter  the  period  of  time,  the  more 
approximate  the  relationship  (eq  I)  because  of  heat  temporarily  stored  within  the  construction. 

Mathematical  understanding  of  eq  1  may  be  obtained  by  paraphrasing  and  simplifying  a  deriva¬ 
tion  by  Sonderegger  (1980).  We  can  characterize  a  construction  thermally  with  a  frequency- 
dependent  matrix  such  that  the  temperature  and  heat  flux  on  an  inside  surface  can  be  described  as  a 
function  of  the  outside  temperature  and  flux: 


3 

Cin(w)  < 


A(u>)  B(cj) 


C(u)  D(  u) 


Ku,  («*>) 

0OUS  (w) 


where  tin  and  tout  are  Fourier  transforms  of  indoor  and  outdoor  temperatures  and  @iB  and 
Qout  are  indoor  and  outdoor  heat  fluxes,  respectively.  A(w),  fl(w),  C(w),  and  D(u)  would  result 
from  matrix  multiplication  of  each  homogeneous  layer  represented  by  matrix  elements  of  the 
form: 

A(cj)  *  cosh  (</•/■( oj))  (A2a) 


^  '  (***M) 

C(w)  *  **F(w)  sinh  (d*F{u>)) 


D( «)  =  cosh  (d*F(ui)) 

F{u)  *  •jlZJJa 


where: 


w  »  angular  frequency  (rad/hr) 

d  •  thicknesa  of  layer  (m) 

i » the  imaginary  unit 

a  ■  diffus'd ty  of  the  layer  (m2  /hr) 

k  >  thermal  conductivity  of  the  layer  (W/(m  k)) 

If  we  rearrange  eq  A2  in  terms  of  we  obtain: 
dm  •  {C-DAfB)tomi  ♦  (D/B)  tm 
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which,  for  the  purely  resistive  indoor  air  layer  has  values  of  A  =  1 ,  B  =  d/k  (=  R-value),  C  =  0,  and 
D  =  l .  This  gives  the  admittances  of  the  wall: 


K,  (oj)  =  D(u>)/B(ui) 

(A4) 

Ym(co)=  1  /fl(w) 

(A5) 

and  a  recharacterization  of  eq  A3: 

Qin  =  r,  My«)  -  rm(«)rou,(w)  (A6) 

In  his  paper,  Sonderegger  treats  fluctuations  in  diurnal  outside  temperatures  and  hourly  cycling 
of  the  heating  plant  as  sine  waves.  Since  the  effect  of  the  heating  plant  cycles  on  the  accuracy  of 
eq  1  is  negligible,  it  is  omitted  in  the  following  discussion. 

Assume  that,  with  sinusoidal  variation  in  outside  temperature  with  a  frequency  of  wout,  an 
amplitude  of  /lou,,  and  a  phase  angle  at  t  =  0  of  7out,  we  can  re-express  outside  temperature  as: 


7-„u.(')=rou,  +  i4ou»  *sin  <Wout  *''Tou«)  (A?) 

Then,  where  rout  is  the  average  value  of  Tout(t),  with  a  constant  Tjn  and  eqs  A6  and  A7,  the 
indoor  heat  flux  variation  can  be  expressed  as: 

<2ln(r)  =  Y,  (0)Tin  -  ym(0)Fou,  -  |rm(«)Mout  («„»,  *7  -  r„ut  +  yj  (as> 

where  ym  is  the  phase  angle  of  Ym  (wout).  For  wCnt  =  0, 

.  (0)  =  Ym(0)-Uo  (A9) 

where  UQ  is  the  steady-state  {/-value  of  the  construction.  Simplifying  further,  we  can  define  t  -  0 
such  that  the  outdoor  temperature  sine  curve  crosses  zero  at  that  time,  so  >ou,  =  0.  This  gives 

Qin(t)  -  U0*AT  -  |Fm(wouf)|/lout  *sin  (wout*f  +  yj  (A10) 

where  AT  =  Tin  -  Tout. 

Equations  A7  and  A10  can  be  rewritten  as: 

Tin  ~  W 0=  AT  11  -Mou|  *sin (w01ll  •DVAT )  (All) 

(Kt) *  U0  *AT  (1  -(|yin(wout)Mou,  *sin  (wout  *t  ♦  ym)/U0*Af)  (M2) 


Now  we  are  ready  to  look  at  the  source  of  eq  1  by  simply  integrating  eqs  A1 1  and  A12  and 
obtaining  their  ratio: 

•i 

(J’a-W'M  _  J»  0 

7 -  "*•  ~T' -  {> 

/  m))dt  f  0  -D)dt 

o  o 

where  R0  ■  l/C/0 

Ar«A0.,*sin(Wo,t*r)  _ 

Da  oat  •sin  7.  VAf. 
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It  is  clear  that  N  and  D  will  always  be  bounded  by  their  coefficients,  whereas  the  rest  of  the 
integral  will  grow  with  time  and  dominate  both  numerator  and  denominator  to  approximately 
the  ratio  of  one.  Therefore,  the  longer  the  duration  of  measurement  t\  the  more  closely  R  (t') 
approximates  R0.  This  is  the  basis  for  demonstrating  the  validity  of  eq  I ,  which  is  simply  a  sum 
approximating  eq  A1 3. 

Actual  diurnal  temperature  cycles  are  not  sinusoidal  or  even  of  regular  amplitude.  However,  if 
one  were  to  decompose  a  segment  of  weather  history  into  its  component  cycles  in  a  Fourier  series, 
each  component  could  be  analyzed  with  eq  A13.  The  real  impact  on  the  calculation  of  eq  A13 
from  imperfect  cycles  would  be  the  length  of  time  required  for  N  and  D  to  become  small  com¬ 
pared  with  the  rest  of  the  numerator  and  denominator  when  integrated,  and  for  R  (r)  to  approxi¬ 
mate  R0  satisfactorily.  The  steadily  diminishing  oscillations  of  R  ( t )  about  a  constant  value  should 
signal  when  a  satisfactory  history  of  data  has  been  obtained. 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 
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